The Blume-Emery-Griffiths Spin Glass and Inverted Tricritical Points 



V. Ongun Ozgelik 1 and A. Nihat Berker 2,3 ' 4 
1 Department of Physics, Istanbul Technical University, Maslak 344-69, Istanbul, Turkey, 
of Sciences and Arts, Kog University, Sariyer 34450, Istanbul, Turkey, 



Department of Physics, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, U.S.A., and 
4 Feza Giirsey Research Institute, TUBITAK - Bosphorus University, Cengelkdy 34684, Istanbul, Turkey 

The Blume-Emery-Griffiths spin glass is studied by renormalization-group theory in d — 3. The 
boundary between the ferromagnetic and paramagnetic phases has first-order and two types of 
second-order segments. This topology includes an inverted tricritical point, first-order transitions 
replacing second-order transitions as temperature is lowered. The phase diagrams show disconnected 
spin-glass regions, spin-glass and paramagnetic reentrances, and complete reentrance, where the 
spin-glass phase replaces the ferromagnet as temperature is lowered for all chemical potentials. 
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The Blume-Emery-Griffiths (BEG) model [H, Q is the 
simplest system for the study of the various meetings 
of first- and second-order phase boundaries between or- 
dered and disordered phases, in a plethora of phase dia- 
gram topologies |3j]. In these diagrams, the second-order 
phase transitions are dominated by thermal fluctuations 
and occur at high temperatures. The first-order phase 
transitions evolve, to finite temperatures, from zero- 
temperature ground-state energy crossings and occur at 
low temperatures. In a well-known phase diagram topol- 
ogy, a tricritical point separates the high-temperature 
second-order boundary and the low-temperature first- 
order boundary. The BEG model has been used to de- 




scribe 3 He- 4 He mixtures [lj, solid- liquid-gas systems 
multicomponent fluid and liquid-crystal mixtures Q , mi- 
croemulsions @, semiconductor alloys @, [H, and elec- 
tronic conduction systems @. 

The inclusion of frozen disorder (quenched random- 
ness) to these systems should yield new phase diagrams, 
as is indeed seen in our current work. We find that 
a temperature sequence of transitions that is reverse 
to the above can occur with the inclusion of quenched 
randomness. Thus, an inverted tricritical point is ob- 
tained, separating a high-temperature first-order bound- 
ary and a low-temperature second-order boundary. Since 
the BEG model is the generic model for tricriticality, 
we believe that the quenched randomness effect of in- 
verted tricriticality should be quite generally a pplic able. 
With frustrated quenched randomness [13, EH, HH EH , 
the spin-glass phase appears within the Blume-Emery- 
Griffiths global phase diagram (Fig.l). Thus, a new 
spin-glass phase diagram topology is found, in which dis- 
connected spin-glass regions occur close to the ferromag- 
netic and antifcrromagnetic phases, but are separated by 
a paramagnetic gap. No spin-glass phase diagram has to 
our knowledge to-date yielded disconnected spin-glass re- 
gions. However, if the global phase diagram of a physical 
realization of the BEG spin glass is fully explored, these 
disconnected regions should be found. 

We have studied, in spatial dimension d = 3, the model 
with Hamiltonian 
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FIG. 1: (Color on-line) Our calculated global phase diagram 
for K = 0. The ferromagnetic phase is bounded by a first- 
order surface close to p = 0, which recedes along the full line 
on the surface from a new second-order transition induced by 
randomness and controlled by a strong-coupling fixed distri- 
bution. At the dashed line, an ordinary second-order transi- 
tion takes over. The transitions from the spin-glass phase, to 
the paramagnetic or ferromagnetic phase, are second order. 
The system being symmetric about p = 0.5, the antiferromag- 
netic sector is not shown. 



where s,* = 0, ±1 at each site i of the lattice and < ij > 
indicates summation over nearest-neighbor pairs of sites. 
The spin-glass type of quenched randomness is given by 
each local Jy being ferromagnetic with the value + J with 
probability 1 — p and antifcrromagnetic with the value 
—J with probability p. Under the scale change induced 
by renormalization-group transformation, all renormal- 
ized interactions become quenched random and the more 
general Hamiltonian 
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FIG. 2: The d = 3 hierarchical lattice for which our calcula- 
tion is exact is constructed by the repeated imbedding of the 
graph as shown in this figure. This hierarchical lattice gives 
very accurate results for the critical temperatures of the d 
=3 isotropic and anisotropic Ising models on the cubic lattice 

M. 



-PH = ]T [JijSiSj+Kyitfj-Aijtf+ty-^tf-g})], 

<tj> 

(2) 

has to be considered. The renormalization-group flows 
are in terms of the joint quenched probability distribution 
P ( Jjj , Kg , Ajj , A j,- ) , which is renormalized through the 
convolution PS 



P'(K r ) = J[]ldK ij P(1^)]8(K' i , j ,-B.({K ij })), (3) 

where primes refer to the renormalized system, Ky = 
(Jij,K{j,Aij,Alj), and R(Ky ■) is the local recursion re- 
lation through which 108 unrenormalized local interac- 
tions in {Kjj} determine 4 renormalized local interac- 
tions in K^/y. The local recursion relation R(Kjj) is 
effected by a mixed Migdal-Kadanoff procedure [15| with 
d = 3 and length rescaling factor b = 3 necessary for 
the equal a priori treatment of ferromagnetism and an- 
tiferromagnetism. Thus, our treatment is approximate 
for the cubic lattice and exact for the hierarchical lattice 
P,[ll[13,[ll,[li,[23,[2ll,[ll[lip shown in Fig.2. This 
hierarchical lattice is known to give very accurate results 
for the critical temperatures of the d =3 isotropic and 
anisotropic Ising models on the cubic lattice [HI]. Fur- 
thermore, in general, exact calculations for hierarchical 
lattices have been seen to constitute very good approxi- 
mations for cubic lattices. 0, I2I I25I I2H \2% 

The probability distribution P(Jy, -Ky, Ay, Ay) is 
represented by histograms lodged on a four-dimensional 
interaction space (Jij,Kij, Ay, Ay). Eq.(3) is effected 
by 8 pairwise convolutions, which are either bond-moving 
or decimation in the appropriate sequence, between inter- 
mediate distributions. The number of histograms rapidly 
grows from the starting two described after Eq.(l). Thus, 
for calculational purposes, before each pairwise convo- 
lution, the histograms are combined by using a bin- 
ning procedure, so that our results are obtained by the 
renormalization-group flows of 22,500 histograms. 

Tricritical phase diagram cross-sections of the purely 
ferromagnetic system for different K/J values are shown 
in Fig. 3. These are standard tricritical phase diagrams, 
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FIG. 3: (Color on-line) Tricritical phase diagram cross- 
sections of the purely ferromagnetic system for the Kj J values 
indicated on the figure, shown consecutively from the inner- 
most curve for K/J = 0. First- and second-order transitions 
are shown by dotted and full lines, meeting at a tricritical 
point. In these systems, with no quenched randomness, the 
standard tricritical topology occurs, with the second-order 
boundary at high temperature and the first-order boundary 
at low temperature. 



in the absence of quenched randomness, with the tricrit- 
ical point separating the second-order transitions at high 
temperature and the first-order transitions at low tem- 
perature. The humped boundary, occurring in mean-field 
theory but not in the d = 2 system is thus found to 
occur in the d — 3 system. 

Our calculated global phase diagram for the BEG spin- 
glass system is in Fig.l for K = 0. The ferromagnetic 
phase is bounded by a first-order surface close to p = 0, 
which recedes along the full line on the surface from a 
new second-order transition induced by randomness and 
controlled by a strong-coupling fixed distribution. At the 
dashed line, an ordinary second-order transition takes 
over. The full line is thus a line of random-bond tricrit- 
ical points. The dashed line is a line of special critical 
points around which universality is violated, since the 
second-order phase transitions on each side of this line 
have different critical exponents. [3| These two lines meet 
at the non-random (p = 0) tricitical point. The transi- 
tions from the spin-glass phase, to the paramagnetic or 
ferromagnetic phase, are second order. 

Cross-sections of this global phase diagram for con- 
stant chemical potential A/ J of the non-magnetic state 
are in Fig. 4. The outermost cross-section has A/J = 
—00, meaning no = states, and therefore is the 
phase diagram of the spin-1/2 Ising spin glass, showing 
as temperature is lowered the paramagnet-ferromagnet- 
spinglass reentrance [25L [26[. The annealed vacancies, 
namely the nonmagnetic states Sj = 0, are introduced 
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FIG. 4: (Color on-line) Blume-Emery-Griffiths spin-glass 
phase diagrams: Constant A/J cross-sections of the global 
phase diagram in Fig.l. The outermost cross-section has 
A/J — —oo, meaning no Si = states. The successive cross- 
sections, going inwards from the outermost cross-section, are 
for the successively higher values of A/J indicated on the fig- 
ure. Thus, annealed vacancies Si = are introduced in these 
cross-sections with successively higher values of A/J, mak- 
ing all ordered phases recede. The dotted and full lines are 
respectively first- and second-order phase boundaries. The 
dashed lines are strong-coupling second-order phase bound- 
aries induced by quenched randomness. The inverted tricrit- 
ical topology is seen between the dotted and dashed lines, 
with the first-order transitions occurring at high temperature 
and the second-order transitions occurring at low tempera- 
ture, on each side of the tricritical point. A new spin-glass 
phase diagram topology is obtained for A/J = 0.35, in which 
the spin-glass phase occurs close to the ferromagnetic (and, 
symmetrically, antiferromagnetic, not shown here) phase, but 
yields to the paramagnetic phase as p is increased towards 
0.5. The spin-glass phase disappears at A/J = 0.37. The 
insets, with expanded scales, clearly show the ferromagnetic 
to paramagnetic phase reentrance and the ferromagnetic to 
spin-glass phase reentrance. 



in cross-sections with successively higher values of A/J. 
For A/J greater than the non-random tricritical value of 
A/J = 0.192, first-order transitions between the ferro- 
magnetic and paramagnetic phases are introduced from 
the low randomness side, but are converted to the strong- 
coupling second-order transition at a threshold value of 
randomness p. This constitutes an inverted tricritical 
point, since the phase boundary is converted from first 
order to second order as temperature is lowered, con- 
trary to the ordinary tricritical points (as seen for ex- 
ample in Fig. 3). The above results are consistent with 
the general prediction that, in d = 3, quenched random- 
ness converts first-order boundaries into second order, 
at a threshold amount of randomness. [27j (In d = 2, 
this conversion is predicted to happen with infinitesimal 
quenched randomness. [27l. [28j ) . The randomness thresh- 
old in d — 3 is of course higher for stronger first-order 
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FIG. 5: (Color on-line) Spin-glass Blume-Emery-Griffiths 
phase diagrams: Constant p cross-sections of the global phase 
diagram in Fig.l. The dotted and full lines are respectively 
first- and second-order phase boundaries. The dashed lines 
are strong-coupling second-order phase boundaries induced 
by quenched randomness. The outermost curve in (a) cor- 
responds to the pure Blume-Emery-Griffiths model with no 
quenched randomness [p = 0). The successive cross-sections 
in (a), going inwards from the outermost cross-section, are for 
the successively higher values of p indicated on the figure. As 
spin-glass quenched randomness is introduced with increasing 
values of p, ordered phases and first-order phase transitions 
recede. 



transitions. In the current system, increased frustrating 
quenched randomness has two parallel effects, namely 
driving the phase transition to lower temperature and 
reaching the threshold for the conversion to second order, 
which explains the calculated results of inverted tricrit- 
icality. (With non- frustrating quenched randomness, on 
the other hand, the transition can actually be driven to 
higher temperature, while the conversion to second order 
still happens (29|.) 

As the annealed vacancies Sj = are increased, at 
A/J > 0.34, of the second-order transitions between the 
ferromagnetic and paramagnetic phases, only the strong- 
coupling transition remains. At A/J > 0.42, the strong- 
coupling second-order transition also disappears, leav- 
ing only first-order transitions between the ferromagnetic 
and paramagnetic phases. Also as the annealed vacan- 
cies are increased, all ordered phases recede. In this pro- 
cess, first the spin-glass phase disappears, at A/J = 0.37, 
which is understandable, since it is tenuously ordered due 
to frustration. The new, disconnected spin-glass phase 
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FIG. 6: Projections of the fixed distributions P*(Jij, ify, Ay, Ajj) for: (a) the disorder-induced second-order transitions 
between the ferromagnetic and paramagnetic phases, (b) the first-order transitions between the ferromagnetic and paramagnetic 
phases, (c) the second-order transitions between the ferromagnetic and spin-glass phases, (d) the second-order phase transitions 
between the spin-glass and paramagnetic phases, and (e) the sink fixed distribution for the spin-glass phase. Note that 
(a),(b),(c),(e) are runaways, in the sense that the couplings renormalize to infinity while the distribution retains its shape 
shown here. In the second-order phase transitions between the spin-glass and paramagnetic phases (d), A is a runaway (to 
minus infinity), while the other interactions remain finite. The fixed distributions in this figure are singly unstable, except for 
the sink (e), which is totally stable. 



diagram topology is obtained in this neighborhood, e.g., 
for A/ J = 0.35 as shown in Fig. 4, in which the spin-glass 
phase occurs close to the ferromagnetic (and, symmetri- 
cally, antiferromagnetic, not shown in the figures) phase, 
but yields to the paramagnetic phase as p is increased to- 
wards 0.5 . The disconnected spin-glass phase diagrams 
occur in a very narrow portion of the global phase dia- 
gram. 

The paramagnetic-ferromagnetic-spinglass reen- 
trances, as temperature is lowered, of the Bhunc- 
Emery-Griffiths spin-glass cross-sections fall on the same 
reentrant second-order boundary, as seen in Fig. 4. As 
seen for A/ J = 0.45 and 0.48 in this figure, before disap- 
pearing at A/J = 0.5, the ferromagnetic phase exhibits 
paramagnetic-ferromagnetic-paramagnetic reentrance as 
temperature is lowered. 

Constant p cross-sections of the global phase dia- 
gram in Fig.l are shown in Fig. 5. The outermost curve 
corresponds to the pure Blume-Emery-Griffiths model 
with no quenched randomness (p = 0). As spin-glass 
quenched randomness is introduced with increasing val- 
ues of p, we see that the first-order boundary recedes 
to the strong-coupling second-order boundary, while the 
ordinary second-order boundary also expands. At p — 
0.18, the first-order transition completely disappears. At 
p = 0.241, the spin-glass phase appears below the fer- 
romagnetic phase, reflecting complete reentrance. At 
p = 0.249, the spin-glass phase completely replaces 
the ferromagnetic phase as the ordered phase, which is 



enveloped by second-order transitions only. Thus, for 
0.241 < p < 0.759, the second-order boundary between 
the spin-glass and paramagnetic phases reaches zero tem- 
perature. 

In the results above, the phase diagrams are deter- 
mined by the basins of attraction of the renormalization- 
group sinks, namely the completely stable fixed points 
and fixed distributions: Each basin is a thermodynamic 
phase. The nature of the phase transitions is determined 
by analysis of the unstable fixed points and fixed dis- 
tributions to which the phase diagram points of these 
transitions flow. Fig. 6 shows the unstable fixed distri- 
butions of (a) the quenched randomness-induced second- 
order transitions between the ferromagnetic and param- 
agnetic phases, (b) the first-order transitions between the 
ferromagnetic and paramagnetic phases, (c) the second- 
order transitions between the ferromagnetic and spin- 
glass phases, and (d) the second-order transitions be- 
tween the spin- glass and paramagnetic phases. The (to- 
tally stable) sink fixed distribution of the spin-glass phase 
is also shown, in (e). The eigenvalue exponent of the 
unstable fixed distribution controlling (b) the first-order 
transitions between the ferromagnetic and paramagnetic 
phases is y = 3 = d, as is required for first-order transi- 
tions [30j . The eigenvalue exponents of the other unstable 
fixed distributions, (a),(c),(d), are y < d as is required 
for second-order transitions. 

We thank C. Guven M. Hinczewski, C.N. Kaplan, and 
O.S. Sanyer for comments. This research was supported 



■5 



by the Scientific and Technological Research Council 



[1] M. Blume, V.J. Emery, and R.B. Griffiths, Phys. Rev. A 
4, 1071 (1971). 

[2] A.N. Berker and M. Wortis, Phys. Rev. B 14, 4946 
(1976). 

[3] W. Hoston and A.N. Berker, Phys. Rev. Lett. 67, 1027 
(1991). 

[4] J. Lajzerowicz and J. Sivardiere, Phys. Rev. A 11, 2079 
(1975). 

[5] J. Sivardiere and J. Lajzerowicz, Phys. Rev. A 11, 2090 

(1975); 11, 2101 (1975). 
[6] M. Schick and W.-H. Shih, Phys. Rev. B 34, 1797 (1986). 
[7] K.E. Newman and J.D. Dow, Phys. Rev. B 27, 7495 

(1983). 

[8] T. Tanaka and T. Kawabe, J. Phys. Soc. Jpn. 54, 2194 
(1985). 

[9] S.A. Kivelson, V.J. Emery, and H.Q. Lin, Phys. Rev. B 
42, 6523 (1990). 

[10] H. Nishimori, Statistical Physics of Spin Glasses and In- 
formation Processing (Oxford University Press, 2001). 

[11] J.J. Arenzon, M. Nicodemi, and M. Sellitto, J. Physique 
I 6, 1143 (1996). 

[12] M. Sellitto, M. Nicodemi, and J.J. Arenzon, J. Physique 
I 7, 945 (1997). 

[13] J.J. Arenzon, F. Ricci-Tersenghi, and D. Stariolo, Phys. 
Rev. E 62, 5978 (2000). 

[14] A. Falicov and A.N. Berker, Phys. Rev. Lett. 76, 4380 
(1996). 

[15] A. Erbas,, A. Tuncer, B. Yiicesoy, and A.N. Berker, Phys. 
Rev. E 72, 026129 (2005). 



(TUBITAK) and by the Academy of Sciences of Turkey. 



[16] A.N. Berker and S. Ostlund, J. Phys. C 12, 4961 (1979). 
[17] R.B. Griffiths and M. Kaufman, Phys. Rev. B 26, 5022R 
(1982). 

[18] M. Kaufman and R.B. Griffiths, Phys. Rev. B 30, 244 
(1984). 

[19] M. Hinczewski and A.N. Berker, Phys. Rev. E 73, 066126 

(2006) . 

[20] M. Hinczewski, Phys. Rev E 75, 061104 (2007). 

[21] Z. Zhang, L. Rong, and S. Zhou, Physica A 377, 329 

(2007) . 

[22] H.D. Rozenfeld, S. Havlin, and D. ben-Avraham, New J. 

Phys. 9, 175 (2007). 
[23] E. Khajeh, S.N. Dorogovtsev, and J.F.F. Mendes, Phys. 

Rev. E 75, 041112 (2007). 
[24] C. Monthus and T. Garel, Phys. Rev. B 77, 134416 

(2008) . 

[25] G. Migliorini and A.N. Berker, Phys. Rev. B 57, 426 
(1998). 

[26] F.D. Nobre, Phys. Rev. E 64, 046108 (2001). 

[27] K. Hui and A.N. Berker, Phys. Rev. Lett. 62, 2507 

(1989); erratum 63, 2433 (1989). 
[28] M. Aizenman and J. Wehr, Phys. Rev. Lett. 62, 2503 

(1989); erratum 64, 1311 (1990). 
[29] A. Malakis, A.N. Berker, LA. Hadjiagapiou, and N.G. 

Fytas, Phys. Rev. E 79, 011125 (2009). 
[30] B. Nienhuis and M. Nauenberg, Phys. Rev. Lett. 35, 477 

(1975). 



